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Abstract 

We study a new non-local hierarchy of equations of the isentropic gas dynamics type where the 
pressure is a non-local function of the density. We show that the hierarchy of equations is integrable. 
We construct the two compatible Hamiltonian structures and show that the first structure has three 
distinct Casimirs while the second has one. The existence of Casimirs allows us to extend the flows 
to local ones. We construct an infinite series of commuting local Hamiltonians as well as three 
infinite series (related to the three Casimirs) of non-local charges. We discuss the zero curvature 
formulation of the system where we obtain a simple expression for the non-local conserved charges, 
which also clarifies the existence of the three series from a Lie algebraic point of view. We point out 
that the non-local hierarchy of Hunter-Zheng equations can be obtained from our non-local fiows 
when the dynamical variables are properly constrained. 



1 Introduction: 



The dynamics of a gas are described by the equations P 



Ut + UUx + - Px = 

P 

Pt + {pu)x = 



St + USx = 



(1) 



where u, p, s denote respectively the velocity, density and entropy while P = P{p, s) represents 
the pressure. For an isentropic gas (s constant), these equations have been well studied in the 
literature for the cases where the pressure is a local monomial of the density. For example, for 
P = , 7 7^ 0, 1, Eq. is known as the polytropic gas equations, and for the choice P = — 
the system is called the Chaplygin gas. This class of systems are known to be of hydrodynamic 
type [D and are integrable. Various properties associated with such systems have been derived in 
the past several years IH [3] • 

In this paper, we will study a new class of equations of the gas dynamics type where the pressure 
is a non-local function of the density. In particular, we will study the system for the choice 
P = -\{d-^vf, where V = p (we make this identification to be consistent with the conventional 
choice in the literature), 



and show that the hierarchy of equations associated with this system is integrable. Even though 
the system belongs to the class of gas dynamics equations, we show that it is not of hydrodynamics 
type and, correspondingly, the proof of integrability is different. The system is bi-Hamiltonian 
[HI Hj (possesses two compatible Hamiltonian structures) and the infinite set of (non-local) charges 
can be constructed recursively. The two Hamiltonian structures of the system possess nontrivial 
Casimirs (distinguished functionals) leading to three distinct infinite series of non-local charges and 
allowing us to extend the flows to local ones [HI El- As a result, we can also construct recursively 
an infinite set of local conserved charges associated with the system which are in involution. The 
zero curvature formulation for both the local and the non-local equations of the hierarchy can be 
given in terms of the Lie algebra valued potentials belonging to SL{2) ® U{\) H) and generate 
the non-local charges in a compact form. However, we have not yet succeeded in finding a scalar 





(2) 
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Lax description for the local equations and, consequently, the relation between the local conserved 
charges and a Lax operator remains an open question. Finally, we show that the hierarchy of 
equations associated with reduce to the non-local hierarchy of Hunter-Zheng equations 
when properly constrained, whereas the local equations become trivial (they do not go over to the 
hierarchy of Harry Dym equations ^2] ) • We note that even though we do not yet know of a physical 
system corresponding to equations ©, the nice integrability properties that emerge make it worth 
studying it in its own right. Furthermore, from its relation with the Hunter-Zheng equation, it is 
quite likely that such a system of equations would find application in a physical problem. 

Our paper is organized as follows. In section 2, we show that the system is bi-Hamiltonian. We 
obtain the two Hamiltonian structures associated with the Hamiltonian description of equation @. 
We prove Jacobi identity as well as the compatibility of the two structures using the method of 
prolongation. We comment on the relation between these structures and the known Lie algebras. 
The bi-Hamiltonian nature of the system is sufficient to guarantee the integrability of the system. 
We construct the recursion operator and obtain the first few charges recursively. We show that the 
two Hamiltonian structures have non-trivial Casimirs which allows us to construct three distinct 
series of non-local charges associated with the three diff^erent Casimirs of the first Hamiltonian 
structure. The Casimir of the second Hamiltonian structure, on the other hand, allows us to 
extend the fiows to local ones. In this case, the recursion operator can be inverted and allows us to 
construct the infinite set of local charges associated with the system recursively. By construction 
these charges are in involution. We also describe briefiy the matrix Lax description for the system. 
In section 3, we present the zero curvature description of both the local as well as the non-local 
hierarchies of equations based on the Lie algebra SL{2) ® U{1). We obtain the non-local charges 
associated with the system from the zero curvature which also clarifies from a different perspective 
why there are three distinct series of non-local charges. In section 4, we summarize our results 
as well as present some open problems. In particular, we point out how the non-local hierarchy 
of equations associated with ((2)) reduce to the non-local hierarchy of Hunter-Zheng equations. In 
spite of having a zero curvature formulation of the system, we have not succeeded in obtaining a 
scalar Lax description for the local equations which remains an interesting open problem. 



2 Bi-Hamiltonian Structure: 



In this section, we will show that equation is a bi-Hamiltonian system with an infinite number 
of conserved charges. From the structure of the equations, it is easy to construct the first few 
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conserved charges of the system which take the forms 

Hi = J dxv , 



H2 



dxuv , 
1 



dx 



1 



2 2^ ' 



(3) 



From the structm'e of the conserved charges in we see that equation can be written in the 
Hamiltonian forms as 

/ u\ / 5H^/5u \ / 5H2/5u \ 

\v V SH3/5V J ' V J ' 



where we have defined 



Pi 



-d 
-d 



V2 



d ^ —Ux 
Ux —{vd + dv) 



(4) 



The skew symmetry of these Hamiltonians structures is manifest. The proof of the Jacobi iden- 
tity for these structure as well their compatibility can be shown through the standard method of 
prolongation which we describe briefly. Introducing the matrix uni-vector 



we can construct the two bivectors associated with the two structures Vi and P2 as 



e^,, = dx AViO^ = - j dx9iA92x 



■| = i y dx 6*1 A {d-^Oi) - 2?; 6*2 A 92x - 2ux 6*1 A 62 



Using the prolongation relations for any vector field v, 



2x ) 
— Olx , 

{d-^ei)-uxe2 , 

Ux Oi -V92x - {v02)x , 



(5) 
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it is straightforward to show that the prolongations of the bivectors Gxii and Bpj vanish, 



implying that 2?i and P2 satisfy Jacobi identity. Furthermore, using ©, it also follows that 



This shows that Pi and P2 are compatible, namely, not only are Pi and P2 genuine Hamiltonian 
structures, any arbitrary linear combination of them is as well. As a result, the dynamical equations 
in ^ correspond to a bi-Hamiltonian system and, consequently, are integrable [SlIZ]. 

It is worth making a few remarks about these Hamiltonian structures. We note that the first 
Hamiltonian structure is the standard structure that arises in systems of hydrodynamic type (for 
example, in polytropic gas dynamics 013]). It is the second Hamiltonian structure which normally 
has some interesting connection with Lie algebras. The Lie algebra structure of P2 is not quite 
manifest in the form given in (jlj). However, with a change of basis 



We recognize this to be the Lie algebra of SL{2) ® U{1) without the central charge (for SL{2)). 
This is similar to the algebra in the case of the two boson hierarchy where there is a central charge 
present in the SL{2) algebra JO]. It is also clear from this analysis that we can naturally assign 
the scaling dimensions \v\ = [v] = 2,[u] = [u] + 1 = l,[x] = — 1 (with these assignments [t] = — 1 
for non-local equations and [t] = for local ones as will be clear later), which will be useful later 
in connection with the zero curvature formulation of the system. For completeness, we also note 
that in this new basis, the first Hamiltonian structure takes the form 



We will, however, continue to use, for simplicity, the variables u, v which are conventional in the 
study of such systems. 



U = Ux 




it follows that the second structure can be written as 
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For bi-Hamiltonian system such as in ((^, we can naturally define an associated hierarchy of 
commuting flows through the relation 



u \ / SHn+i/6u \ / 5Hn/5u 



n = 1,2, 



The gradients of the successive Hamiltonians in the hierarchy can be related through the recursion 
operator as 



6Hn+i/Su 
SHn+i/6v 



fit 



SHn/Su 



(6) 



where we have defined 



fi = V2VI 



-1 



Uxd ^ —d 



\ 2v + Vxd ^ —Uxd ^ J 
It is interesting to note that, in this case, we can invert the recursion operator to write 



(7) 



fi 



-1 



-^dzd ^zuxd 



\dzd 



where we have defined 



-d^ — ^d'^zuxd ^Uxzd ^d'^u^zd ^z 



I 1 2\-l/2 

z=[y- ^Ux) 



(8) 



This will be useful later in the construction of the local flows associated with the system. Using 
Eqs. © and O, we can write the relations for the gradients explicitly as 



6H, 



-d-^ux—^ + {2v - d'^vx)- 
ou du dv 

5Hn , 6Hn 



"+1 ^ _9-2!^ + a-V^, n=l,2,.... 
ou ov 



6v 

These can be explicitly integrated to give the infinite set of (non-local) Hamiltonians 
Hi = J dxv , 

H2 = dxuv , 



(9) 



dx 



1^2 1(^-1 )2 

2 2^ ^ 



Ha 



I 



dx 



dx 



-u'^v — uv(d '^v) — -u(d '^v)'^ 
6 2 



— u'^v + 
24 2 



(10) 



The corresponding flows (the first few) have the forms 



1 



{uv)x 



-uvux — -u^Vx + 2v{d ^v) + Vx{d ^w) 



(11) 



We note that, with the dimensionalities of the variables described earlier, all the conserved charges 
in (|lUj) have the same canonical dimension of 1 and by construction they are all in involution. 

To further understand the properties of this hierarchy, let us note that the first Hamiltonian 
structure in ((IJ) has three Casimirs 



dx V 



H 



(1) 



dx Ux ^ Q , 



H? - ] 



dx u 



such that 



/ 6Hi/6u \ I 6H\'> /5u \ I 6Hf^ /5u 



5Hf^ /5v 



(12) 



We remark here parenthetically that the Casimir H^"^ is trivial much like in the case of the Hunter- 
Zheng equation 11 . The existence of Casimirs would normally imply that the series of recursive 
flows cannot be extended to negative values of n. However, in the present case, it is not very hard 
to check that the second Hamiltonian structure P2 also has a Casimir of the form 



H^i = 2 J dx {v-^uiy^'^ , 
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such that 



V2\ =0 

V 6H^i/6v J 



(13) 



As a consequence, the hierarchy of flows in ((TT|) can be extended to negative values of n [S]. We 
note that the Casimir in ()13() is conserved under the flows of (jllj) . 

For negative values of n, the gradients of the Hamiltonians will satisfy a recursion relation 
involving R^^ given in Q, and take the explicit forms 



6Hn _ 1 1 2^-l/2 1 „,2A "1/2 'J-H'n+l 

2^x 



6v 



5u 



6v 



6u 



6u 



5v 



6v 



6Hn r.2^Hn+i 

oux— o — , n = -2, -3, 



6v 6v 

The corresponding conserved charges can now be recursively constructed and have the forms 



(14) 



H^i = 2 J dx {v-\uiy 

7 / 1 2N-1/2 

dxuxx [v - ^U^) 



H. 



1 

12 



--^ j dx {vxx + 2n^2. - UxUxxx) [v -\ ul) ^^'^ 



-\ I dx {VxxUxx - VxUxxx) {v - ^ul) 



dx 



linlx + {v-\ul)J-^^{{ 



64 



V — t;U: 



) 

'/ XX/ 



XJ XX/ 



192 '^xUxx^ {y 2'^x) XXX ^'^xxUxxxx {y 2^x) 



{V - ^Ux) 



l,.2\-7/2 



(15) 
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The dynamical equations following from these (the first few) take the forms 



1 „,2^-l/2■ 



{v-\ ul) 



1 

2 L 



Uxx \y - 2Ux 



l,.2\-3/2 



{V 



1 „,2^-l/2■ 



2 "^x 



Vx {v - \ul) ^'"^ 



(16) 



We note that these flows and the Hamiltonians for negative n are completely local. The Hamil- 
tonians involve increasing number of derivatives of the variables, unlike the conserved charges in 
the polytropic gas where they are pure polynomials of the dynamical variables. Similarly, the local 
dynamical equations become increasingly more nonlinearly dispersive as n becomes more negative. 
Thus, this hierarchy of equations is very different from the usual polytropic gas systems [H 
Nonetheless, all these local charges are in involution (having been constructed from compatible 
Hamiltonian structures). We note here that, with the dimensionalities for the variables described 
earlier, all the local charges in (|15j) have the scaling dimension 0. 

The three Casimirs in (|12|) can be easily checked to be conserved under the flows (|16|) . In fact, 
we can also construct non-local charges from H^^ recursively and they take the forms 



H 



(1) 



dx Ux —>■ 



r(l) 



H'^^' = - I dx (d-^v) 



Hi^^ = I dx d-^(ux{d-\)'^ , 



H 



(1) 



dx d ^ 



._(d-^yf-Uxd-'^Ux{d-\) 



(17) 



(1) 



All these charges have the dimensionality and are in involution. (Note that we can obtain H2 
from the trivial Casimir h[^^ using the recursion relation Q with the prescription (d^^O) = — 1, 
for more details see Ref. jS| and references therein.) Similarly, also leads to the following 

series of non-local charges which are related recursively. 



HP . f 



H. 



(2) 



dxu , 



dx 



-u^ + {d'\) 
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H. 



(2) 



H, 



(2) 



dx 



dx 



.1^4 _ 1(5-2^)2 _ 1^-2(^-1^)2 _ i„(5-2^2) 

24 2^ ' 2 ^ ' 2 ^ ' 



(18) 



The scaling dimensions for this set of charges turn out to be —1. They are ah conserved under the 
local flows of H16|) and are in involution. The meaning of the three series of non-local charges is 
quite clear from the point of view of the existence of three Casimirs. However, in the next section, 
we will see within the context of the zero curvature formulation that the three series are related to 
the fact that the Lie algebra of SL{2) (related to the second Hamiltonian structure of the system) 
has three generators. We remark here that non-local flows can also be derived from the set of 
charges in H17|) and (|18() . but we do not get into that. 

To close this section, we note that a bi-Hamiltonian system of evolution equations, 



Kn[u,v\ 



is known jl3l I14j to have a natural Lax description of the form 

dM 

^ = [M,B], 



where, we can identify 



Here represents the matrix Frechet derivative of Kn , defined by 




— Kn[u + ewi.v + ew2] 



e=0 



In this way, we can obtain a matrix Lax description (M, B are matrix operators) for the non-local 
as well as the local equations. However, such a Lax description is not very useful since it does not 
directly lead to conserved charges. Therefore, we do not give details of this and study the zero 



10 



curvature formulation for this system in the next section. 



3 Zero Curvature: 



To construct the zero curvature for the local as well as the non-local equations, let us recall some of 
the features of our system of equations. We note that the second Hamiltonian structure corresponds 
to the Lie algebra of SL{2)®U{1) so that the zero curvature condition can be based on this algebra. 
Furthermore, the canonical dimensions of the variables are given by [u] = 0, [v] = 2, [x] = — 1, 
[t] = (for local flows). Since the canonical dimension of t is zero, for a zero curvature condition 
of the form 

atAi -a^Ao- [Ao,Ai] = , (19) 

we note that multiplication with the matrix Aq must preserve the canonical dimensions of the 
elements of any matrix and, therefore, will have the unique form 



Ao = , (20) 




where [ ] represents the dimensionality of the matrix element. It follows from l|19j) that the matrix 
Ai will have to have the form (since [d] = 1) 



Ai = . (21) 




Following the procedure in ^U] which describes the general construction of zero curvature based 
on SL{2) (8) ^^'^(l), and recalling the dimensionalities of our matrices in (|20|) and (|21j) . let us choose 



An 



-\B., + ^{-{d~^A) + 2u.,B) -B 

^{-A + \Bv) ~^{d-^A) 



XUr. 



Ai = I ' 1 , (22) 



where A represents a dimensionless spectral parameter and A, B are arbitrary functions of the 
dynamical variables as well as the spectral parameter. The zero curvature condition ()19|) . in this 
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case, leads to the dynamical equations 

ut = \{-{d-^A)+u^B)-B^ , 
vt = X {-UxA + {vd + dv)B) - A,, . 



(23) 



Since the dynamical variables are independent of the spectral parameter, it is clear that the func- 
tions A, B must depend on A for this equation to be meaningful. Let us make a Taylor expansion 
in A of the form (recall that for the local flows n is negative) 



with 



A. 



l„,2\-l/2 



2 xj 



Substituting (jSU and ((221) into we obtain 

A.; \ „ / Aj+i 



Ut 
Vt 



B 



An+l 
Bn+1 



It is clear now that if we identify 



B-i = {v- \ul) 



l,.2\-l/2 



-1,-2 



, . . . , 



Inl + 1 



B, 



6v ' 



(24) 



(25) 



(26) 



(27) 



then (|26|) gives the dynamical equations of the hierarchy (for any n), the two Hamiltonian structures 
of the system as well as the recursion relations between the conserved charges given in ()14[). The 
first local flow for n = — 1, for example, takes the form 



Ut 



Vt 



l„,2\-l/2 



Ux [V - ^Ux) 



which coincides with the first fiow in (fT6|) . 
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For the non-local equations (n positive) let us choose 

-B, + ^{-{d-'A) + 2u,B) -B 



An 



The zero curvature condition ()19l) . in this case, gives 

ut = {-{d-^A) + u^B) - XB^ , 
vt = {-UxA + {vd + dv)B) - XA^ . 

Using a Taylor expansion of the type in H24|l (here n is positive), namely, 

n n 

with 



Ai = 0, Bi = l, 



^i( t'^' )=^2( ), J = l,2,...,n-1, 



we obtain 



Bj+i J \ Bj 

Vt J \ Bn 

Once again, with the identification in (|27|) . it is easy to see that this leads to the dynamical non-local 
equations of the hierarchy. For n = 1, for example, it gives the chiral boson equation of 1)11(1 . 

As we have seen, the zero curvature formulation is quite nice in that it not only gives the 
dynamical equations, but also leads to the two Hamiltonian structures of the system as well as the 
recursion relation for the conserved charges. It can also give the non-local charges (if present) of 
the system [l"51j . For example, let us consider the matrix Ai in ()22|) . Then, it is straightforward to 
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show that as a result of the zero curvature condition H19|) , 



Z = F{e 



dx Ai 



is conserved. Here "P" stands for path ordering of the exponential. This path ordered exponential 
can be expanded as 

/oo roo POO 

dxAi + / dxAi (d-^Ai) - dxAi {d'^ (Ai {d-^Ai))) +■■■ , 
-oo J —oo J ~OD 

and each term of the expansion will be individually conserved. (Actually, the coefficient of each 
independent power of A would be conserved. However, as we will see below, each term in the 
expansion leads to different powers of A for the matrix elements, as a result of which each term 
in the expansion is individually conserved.) We can work out explicitly the first few terms in the 
expansion to see that 

/ POO 2 

' A / dx 11.^ 

A 



dxAi 



A / dxux 

oo 



dx 



\ 



A3 

y — / dxv 



( 







A3 



dxv 



dxAi [d-^Ai) 



I 



X' / dx(u4d-\,)-Ud-'{v-'^ul)) 



4 foo 



^ / dxv{d 



- r dxu.,{d-H) ^ 

J — oo 
\ 2 roo 

/ dxv{d-h) 

2 J-OO 



/ ^2 /"OO POO \ 

' / dx{d~^v) / dxu 

^ J —oo J —oo 

^4 poo ^2 poo 

— dxuv TT dx{d^^v) 

^ J —oo ^ J —oo ' 



dxAi {d~^ (Ai (a-^Ai))) 



dx {d-^Ai) (Ai (9~^Ai)) 
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A3 
2 



dxd' 



dx 



A3 



V 



These are indeed the three series of non-local charges for our system (up to multiplicative constants). 
The origin of the three infinite series of charges can be understood in this approach in the following 
manner. We note that our zero curvature condition is based on the Lie algebra SL{2) iJD ^(1) and 
that the potential Ai belongs to this algebra. However, the Abelian (U{1)) part of the potential 
is not subjected to path ordering and, consequently, docs not contribute to non-local charges. On 
the other hand, SL{2) has three generators and the charges projected along any of the generators 
must be conserved. Consequently, the system possesses three infinite series of non-local charges. 



4 Conclusion: 



In this paper, we have studied a new system of gas like equations where the pressure is a non-local 
function of the density and have shown that the hierarchy of equations is integrable. We have shown 
that the system possesses two compatible Hamiltonian structures and have constructed an infinite 
number of (non-local) conserved charges recursively. The second Hamiltonian structure of this sys- 
tem corresponds to the centerless SL{2) 'S>U{1) algebra. We have shown that the first Hamiltonian 
structure of the system possesses three nontrivial Casimirs while the second Hamiltonian structure 
has one. This allows us to extend the non-local flows into local ones corresponding to negative flows 
of the hierarchy. We have constructed an infinite series of commuting local charges associated with 
the system recursively. This local system of equations possesses three infinite series of non-local 
charges and we have constructed them from the three Casimirs of the first Hamiltonian structure 
recursively. We have given the zero curvature formulation for the system of local as well as non-local 
equations based on the Lie algebra SL{2) (g) [/(I). This brings out naturally the two Hamiltonian 
structures of the system as well as the recursion relations between the conserved charges. The zero 
curvature formulation also leads to the three infinite series of conserved non-local charges of the 
system in a simple manner and relates their existence to the fact that the Lie algebra of SL{2) has 
three generators. 

There are some other interesting features of this system that we would like to discuss briefly. 
We note that the non-local hierarchy of flows associated with the Hunter-Zheng equation can be 
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obtained from when the dynamical variables are constrained as 

u={d''w), v = ^{d-'wf = lul. (28) 

Under this reduction, the three series of non-local charges go over to the ones associated with the 
Hunter-Zheng equation. However, under this reduction, the local flows of the hierarchy become 
trivial and do not go over to the Harry Dym equation. 

Another interesting issue that remains open is the construction of a scalar Lax representation 
for the system of local equations. (We do not expect to find a scalar Lax representation for the 
non-local flows.) It is well known that given the zero curvature formulation of a system, one can 
easily go over to a scalar Lax description and vice versa I17j. In this case, however, in spite 
of the fact that we have a zero curvature formulation, we have not been able to find a scalar Lax 
description for this system. We have discussed the standard matrix Lax description, but finding a 
scalar Lax operator remains an open question that deserves further study. 
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